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Abstract. QCD thermodynamics is investigated by means of the Polyakov-loop-extended Nambu Jona- 
Lasinio (PNJL) model, in which quarks couple simultaneously to the chiral condensate and to a background 
temporal gauge field representing Polyakov loop dynamics. The behaviour of the Polyakov loop as a function 
of temperature is obtained by minimizing the thermodynamic potential of the system. A Taylor series 
expansion of the pressure is performed. Pressure difference and quark number density are then evaluated 
up to sixth order in quark chemical potential, and compared to the corresponding lattice data. The validity 
of the Taylor expansion is discussed within our model, through a comparison between the full results and 
the truncated ones. 

PACS. 12. 38. Aw General properties of QCD (dynamics, confinement, etc.) - 12.38.Mh Quark-gluon 
plasma 



1 Introduction 

QCD thermodynamics has been the subject of intense 
investigations in recent years. Thanks to lattice simula- 
tions, the equation of state of strongly interacting matter 
is now at hand as a function of temperature T and in 
a limited range of quark chemical potential /i. Improved 
multi-parameter re- weighting techniques Taylor se- 

ries expansion methods |5U3U5| and analytic continuation 
from imaginary chemical potential |g||7||§l|y] provide lattice 
data for the pressure, entropy density, quark density and 
selected susceptibilities. In response to these lattice sim- 
ulations, many phenomenological models have been pro- 

to give an interpretation of the available lattice data in 
terms of effective degrees of freedom, and to explore those 
regions of the phase diagram that cannot be reached on 
the lattice yet. 

In this context, encouraged by the successful descrip- 
tion of N c = 2 QCD obtained in the Nambu Jona-Lasinio 
model we have investigated recently full QCD ther- 
modynamics at zero and finite quark chemical potential in 
the framework of a Polyakov-loop extended Nambu Jona- 
Lasinio (PNJL) model [25l2fil27128l29l3()l31| in which quar 
develop quasiparticle masses by propagating in the chiral 
condensate, while they couple at the same time to a homo- 
geneous background (temporal) gauge field representing 
Polyakov loop dynamics. 

The "classic" NJL model incorporates the chiral sym- 
metry of two-flavour QCD and its spontaneous break- 
down at T < T c . Gluonic degrees of freedom are "inte- 



grated out" and replaced by a local four-point interaction 
of quark colour currents. Subsequent Fierz transforma- 
tions project this interaction into various quark-antiquark 
and diquark channels. The colour singlet qq modes of low- 
est mass are identified with the lightest mesons. Pions 
properly emerge as Goldstone bosons at T < T c . However, 
the local SU (N c ) gauge invariance of QCD is now replaced 
by a global SU(N C ) symmetry in the NJL model, so that 
the confinement property is lost. Consequently, standard 
NJL-type models are bound to fail in attempts to describe 
N c = 3 thermodynamics around T c (and beyond) for non- 
zero quark chemical potential \x. 

The deconfinement phase transition is well defined in 
the heavy-quark limit where the Polyakov loop serves as 
an order parameter. In the presence of dynamical quarks 
the Z(2>) center symmetry of the SU{3) gauge group is 
explicitly broken. No rigorous order parameter is estab- 
lished for the deconfinement transition in this case |32|. 
but the Polyakov loop still serves as an indicator of a rapid 
crossover towards deconfinement. 

In the present paper we explore the capability of an 
updated PNJL model to describe essentials of QCD ther- 
modynamics around and above T c . We use a modified 
ksrersion of our previous PNJL model |3l)U-il| . with an im- 
proved effective potential for the Polyakov loop field. One 
of our main tasks here is to check the validity of the Tay- 
lor expansion approach at finite chemical potential fi. The 
presently available lattice simulations that we want to 
compare our results with, provide the equation of state 
of two-flavor QCD at finite quark chemical potential by 
means of a Taylor expansion around fi = up to sixth 
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order in [i. Our PNJL model enables us to produce both 
the full result at finite [i and the Taylor expansion around 
(j, = 0. In this way we can draw conclusions about the 
range of applicability of the Taylor expansion and its con- 
vergence properties. 



2 The model 

The Lagrangian of the Polyakov-Loop-extended Nambu 
Jona-Lasinio (PNJL) model is written in the following 
way: 



C = $ (i^D» - m ) i> + j 
-W(«P,**,T), 



(i) 



where tp = (ip U7 ipd) is the quark field, too — diag(mo, Too) 
is the isospin-symmetric bare quark mass matrix and we 
define 



D» = d" - iA" with A" =5% A , 



(2) 

with Ao = —iA^\ X a are the eight Gell-Mann matrices. 
A local, chirally symmetric scalar-pseudoscalar four-point 
interaction of the quark fields is introduced with an ef- 
fective coupling strength G. The potential W(0,0*,T) in 
the Lagrangian governs the dynamics of the Polyakov 
loop = Tr c (L)/3 and its conjugate 0* = Tr c (L+)/3. 
The matrix L is written in terms of the temporal gauge 
fields: 



L 



V exp 




exp 



iA A 
~T~ 



(3) 



In a convenient gauge (the so-called Polyakov-gauge) the 
Polyakov loop matrix can be given a diagonal representa- 
tion |27|. Note that in the chiral limit (too — > 0), this 
Lagrangian is invariant under the chiral flavour group, 
SU(2)l x SU{2)r, just like the original QCD Lagrangian. 

The effective potential U (0, 0*,T) must obey the fol- 
lowing general features: it must satisfy the Z(3) center 
symmetry, like the pure gauge QCD Lagrangian; besides, 
in accordance with lattice predictions for the behaviour of 
the Polyakov loop, hi must have an absolute minimum at 
0=0 at small temperatures, while above the critical tem- 
perature (To ~ 270 MeV in pure gauge QCD) the mini- 
mum is shifted to a finite value of 0. In the limit T — > oo 
we have — > 1. In our previous works |.'itJII.'il| we have 
chosen the simplest possible form for 14, namely a poly- 
nomial in 0, 0*. An improved effective potential, follow- 
ing Ref. [57J, replaces the fourth order term by log[J(<P)], 
where .7(0) is the Jacobi determinant which results from 
integrating out six non-diagonal group elements while keep- 
ing the two diagonal ones to represent L. We therefore use 
the following ansatz for hi: 



W(0,0*,T) 



1 



-b 2 (T)0*0 



(4) 



+b 4 (T) log 1 - 60*0 + 4 (0* 3 + <Z> 3 ) - 3 (0*0) 



4 - 



-o 
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Fig. 1. Scaled pressure, entropy density and energy density 
as functions of the temperature in the pure gauge sector, com- 
pared to the corresponding lattice data taken from Ref. [55]. 





ai 


a 2 




3.51 


-2.47 


15.22 


-1.75 



Table 1. Parameter set used in this work for the Polyakov 
loop potential J1]|HJ|. 



ith 



b 2 (T) = a +ai 



'f)+« 2 (f) , h(T) = b 4 (^ 

'(5) 

Through its logarithmic divergence as 0, 0* — > 1, this 
ansatz guarantees that the Polyakov loop will be auto- 
matically constrained to be always smaller than 1, reach- 
ing this asymptotic value as T —* oo. 

Following the procedure of [301, a precision fit of the 
parameters ai, bi is performed to reproduce lattice data 
for pure gauge QCD thermodynamics and the behaviour 
of the Polyakov loop as a function of temperature. The 
results of this combined fit are shown in Figs. ^ and 
(dotted line). The corresponding parameters are listed in 
Table ^ The critical temperature To for deconfinement 
in the pure gauge sector is fixed at 270 MeV, in agree- 
ment with lattice results. The NJL part of the model in- 
volves some free parameters, which we take from Ref. |3()| . 
where they were fixed to reproduce some physical quan- 
tities in the hadronic sector. These parameters, and the 
corresponding physical quantities, are listed in Table El 

Before passing to the actual calculations, we summa- 
rize basic assumptions behind eq. and comment on 
limitations to be kept in mind. In fact the PNJL model JTJ 
is quite schematic in several respects. It reduces gluon dy- 
namics to a) chiral point couplings between quarks, and b) 
a simple static background field representing the Polyakov 
loop. This picture cannot be expected to work beyond a 
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Fig. 2. Dotted line: PNJL fit of Polyakov loop as a function 
of temperature in the pure gauge sector, compared to corre- 
sponding lattice results (empty symbols) taken from Ref. |34| . 
Continuous line: PNJL model prediction of the Polyakov loop 
behaviour as a function of the temperature, in the presence of 
dynamical quarks. The corresponding lattice data (full sym- 
bols) are taken from Ref. [85) . 

limited range of temperatures. At large T, transverse glu- 
ons are known to be thermodynamically active degrees of 
freedom, but they are ignored in the PNJL model. To what 
extent this model can reproduce lattice QCD thermody- 
namics is nonetheless a relevant question. We can assume 
that its range of applicability is, roughly, T < (2 — 3)T C , 
based on the conclusion drawn in Ref. |36| that transverse 
gluons start to contribute significantly for T > 2.5 T c . 



3 Results 



A [GeV] G[GeV~ 2 ] m [MeV] 

0.651 10.08 5.5 

|(^V u )| 1/3 [MeV] /x[MeV] m w [MeV] 
251 92J3 139.3 

Table 2. Parameter set used in this work for the NJL model 
part of the effective Lagrangian JJ), and the resulting physi- 
cal quantities. For these values of the parameters we obtain a 
constituent quark mass m =325 MeV. 



ues. In the present work we stay at the mean-field limit 
which consistently implies <P = <P* . 

Minimizing L?(T, /i, a, <P) determines the chiral con- 
densate a and the Polyakov loop <P as functions of temper- 
ature and quark chemical potential. Fig.[3shows the tem- 
perature dependence of the Polyakov loop at fi = for two 
situations: first for the pure-gauge case (dotted line), and 
secondly for the case including quarks (full line). While 
the pure-gauge lattice results have been used to fix the 
effective potential for the Polyakov loop filed as explained 
previously, the case with inclusion of quarks is a prediction 
of the model without any additional tuning of parameters. 
The resulting agreement with corresponding lattice data 
incorporating dynamical quarks j35j is striking. 

Notice that, while the Polyakov loop shows a discon- 
tinuity at the critical temperature in the pure gauge sys- 
tem indicating a first order deconfinement phase transi- 
tion, this turns into a smooth crossover when quarks are 
introduced. At the same time the critical temperature is 
reduced from To = 270 MeV in the pure gauge case to 
around T c — 215 MeV in the presence of quarks (not evi- 
dent from Fig.fJlsince the horizontal axis is scaled by T c ). 



3.1 PNJL model at mean-field level 



3.2 Comparison with lattice results 



After performing a bosonization of the PNJL Lagrangian 
and introducing the auxiliary fields a and ir the thermo- 
dynamic potential results as follows: 



r2(T,/ 1) ff,#,#*)=W(#,#*,T) 



(6) 



+ 2G~ 2N * 



d 3 P 

(27T) 3 



{3E P 6 (A 2 - p 2 



-Tin 
-Tin 



1 + 3<p e - (Ep - fJ - )/T + ■Z$* e - 2{E v-^ )/T +e- 3( - E p-v)/ T 
l + 3<p* c -( E p+t*)/ T + 3$ c - 2 ( E P+v)/ T +c -HE P +n)/T | 



where the quark quasiparticle energy is E p = \/p 2 + m 2 
and the dynamical (constituent) quark mass is the same as 
in the standard NJL model: m = m — (u) = rn — Gfyi/j). 

In general, the fields and <P* (and their thermal ex- 
pectation values) are different at non-zero quark chemical 
potental jSZ] ■ We demonstrate in a forthcoming paper [SHI 
that <P ^ <P* at fi 7^ is primarily a consequence of quan- 
tum fluctuations of the fields around their mean field val- 



In our previous work |3(J| we have compared predictions 
of the PNJL model concerning QCD thermodynamics, at 
zero and finite fJL, with corresponding lattice data. At non- 
zero chemical potential the lattice "data" are deduced 
from a Taylor expansion of the thermodynamic quantities 
in powers of fi/T around /i = 0. For this reason, following 
also Ref. [35], we perform the same expansion in the PNJL 
model, and compare the lattice expansion coefficients to 
those calculated in our model. We have: 



Cn (T) 



f>m(£)" ; 

71 = 

1 9(p(T,/i)/T 4 ) 



n\ d^/T) n 



(7) 



Our results for the coefficients ci and C4 are shown in 
Fig. [3] Once a limited set of input parameters is fitted to 
Lattice QCD in the pure gauge sector and to pion prop- 
erties in the hadronic sector, the PNJL model evidently 
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Fig. 3. Expansion coefficients C2 and C4: the PNJL model results are compared with the corresponding lattice results taken 
from Ref. 0. 



provides a very good description of the lattice data for 
the expansion coefficients of the Taylor series. The results 
for cq (not shown here) are also in very good agreement 
with the lattice data and will be found in a forthcoming 
paper |4*0| . 

From these coefficients the Taylor-expanded pressure J7J) 
is reconstructed and the quark number density follows as: 



T3 



2 C2 (f) + 4c 4 (^) 3 + 6 C6 (^) ; 



(8) 



Fig. 0] displays results for the pressure difference, Ap = 
p(T,p) —p(T,fi = 0), and for the quark number density, 
both expanded around p. = up to sixth order in quark 
chemical potential, consistently with the corresponding 
lattice data. The full computation of the pressure and 
quark number density at finite chemical potential cannot 
be performed on the lattice, but this is of course possible in 
our model. We can thus investigate convergence properties 
of the Taylor series and study discrepancies between the 
full and truncated results. This study is shown in Fig. |SJ 
The continuous lines (full results) are compared to second 
order (dashed) and fourth order (dotted) expansions in £. 
Evidently the series converges rapidly for small chemical 
potentials: the second order differs by less than 5% from 
the full result, and the fourth order basically coincides 
with the solid curve. It is nevertheless also evident that 
when p/T c becomes larger, of order 1, the Taylor expan- 
sion develops problems around T c . In this region in fact, 
the expansion coefficient C4 is rather large, and for this 
reason the fourth order contribution overshoots the full 
one around T c . When T becomes larger, the coefficient C4 
drops and for this reason the agreement between the full 
result and the truncated one is again very good. 



tures: spontaneous chiral symmetry breaking and confine- 
ment. The respective order parameters (the chiral quark 
condensate and the Polyakov loop) are given the meaning 
of collective degrees of freedom. Quarks couple to these 
collective fields according to the symmetry rules dictated 
by QCD itself. Once a limited set of input parameters 
is fitted to Lattice QCD in the pure gauge sector and 
to pion properties in the hadron sector, the quark-gluon 
thermodynamics above T c up to about twice the critical 
temperature is well reproduced. In particular, the coeffi- 
cients of the Taylor expansion in p/T are in remarkably 
good agreement with the ones coming from lattice sim- 
ulations. The PNJL model correctly describes the step 
from the first-order deconfinement transition observed in 
pure-gauge Lattice QCD to the crossover transition when 
Nj = 2 light quark flavours are added. The predicted 
behaviour of the Polyakov loop as a function of temper- 
ature in the presence of two quark flavours is in surpris- 
ingly good agreement with the corresponding lattice data: 
a highly non-trivial result. 

The Taylor-expanded pressure difference and quark 
number density of the system has been calculated up to 
sixth order in quark chemical potential. The agreement 
with corresponding lattice data is again very good. Fi- 
nally, a comparison between the full result for these two 
quantities and the truncated one has been performed. We 
observe fast convergence of the expansion in powers of p/T 
at small chemical potential. At larger chemical potentials, 
comparable to T c , limitations in the applicability of the 
Taylor expansion become apparent. 
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